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Solutions for a class of wave equations with effective potentials are obtained by a method of 
a Laplace-transform. Quasinormal modes appear naturally in the solutions only in a spatially 
truncated form; their coefficients are uniquely determined by the initial data and are constant only 
in some region of spacetime — in contrast to normal modes. This solves the problem of divergence 
of the usual expansion into spatially unbounded quasinormal modes and a contradiction with the 
causal propagation of signals. It also partially answers the question about the region of validity of 
the expansion. Results of numerical simulations are presented. They fully support the theoretical 
predictions. 



I. INTRODUCTION 



We study a class of wave equations [23( in the form 

il>(t,x)=0, (1) 



92 d2 \r I , 



describing propagation of a real wave function, ip : M + x 
K — * K or C, in presence of a real static potential V(x) 
and satisfying initial conditions 



iP(0,x) = f(x), ij>(Q,x)=g(x). 



(2) 



Equations of this type appear, among others, in the ana- 
lysis of linear perturbations around static solutions of 
nonlinear wave equations 0, • I n particular, they de- 
scribe propagation of scalar, electromagnetic and grav- 
itational perturbations of black holes which are static 
solutions of the Einstein equations (Regge- Wheeler eq. 
Q, Zerilli eq. Q). By a suitable change of variables 
and separation of the angular dependence, the original 
wave equations on a curved geometry can be brought to 
the above form ftfr. which contains only second partial 
derivatives and an effective potential V ||. The range 
of the new variable x is determined by the coordinate 
change and usually extends from — oo (corresponds to 
the black hole horizon) to +oo (corresponds to r — > oo). 
The potential contains all necessary information which is 
needed to describe scattering of the (partial) waves on a 
given geometry. 

Study of the linear perturbations of black holes has led 
to the observation that solutions can be expressed as a 
sum over a countable set of modes, called quasinormal 
modes (to be distinguished from normal modes). This 
representation is universal in the sense that the set of 
modes depends only on the form of the potential and not 
on the initial data. Moreover, the modes possess differ- 
ent damping factors and hence only a few least damped 
modes are observed, what makes the picture simplier and 
even more universal. The parameters of a black hole 
are encoded in the set of modes and may be read off by 
future astrophysical observations of gravitational radia- 
tion (LIGO or VIRGO) coming from the vicinity of black 
holes. 



In the following we will leave the interpretation of 
waves travelling on a black-hole background and study 
the mathematical properties of the equation (JTJ, since it 
has wider applications in various areas of mathematical 
physics. 

Originally, quasinormal modes (QNMs) were defined 
as solutions of an eigenvalue problem resulting from time 
separation in the wave equation l(T[l. after imposition 
of (artificial) boundary conditions of an outgoing wave 
at both "ends": x — ► ±oo _6] . It has been shown for 
several potentials that there exists an infinite countable 
set of such modes, what made them a good candidate 
for representation of solutions. Later, Sun and Price 
and Nollert and Schmidt Q used a Laplace trans- 
form approach and a Green's function technique to show 
how a sum over quasinormal modes appears in the solu- 
tion. The boundary conditions of outgoing waves, which 
seemed not to be uniquely defined (for Re(s) < 0), have 
been replaced by a more familiar one of square inte- 
grability (for Re(s) > 0) and analytic continuation (to 
Re(s) < 0), and have been reproduced in the calcula- 
tion in a uniquely defined way. Although the technique 
has been successfully applied for the general definition of 
QNMs and for indication how they appear in the solu- 
tion, it has never been studied (to our best knowledge) 
more rigorously in a general setting, in particular, to clar- 
ify in which region of spacetime the quasinormal modes 
describe an exact solution and what is the form of the 
solution outside the region. There exists only literature 
concerning special potentials H [H EJ 

Moreover, 

dealing with QNMs leads to several mathematical prob- 
lems, which require better understanding of the the role 
they play in the solution. 

First of them is that the profiles of QNMs are spatially 
unbounded. It makes every series built of them divergent 
at big distances. Although, the time dependence of the 
quasinormal modes is different than by normal modes - 
they oscillate and decay, what makes the divergence more 
harmless for late times. Yet the problem of early and in- 
termediate times remains open. It causes also difficulties 
in the definition of excitation coefficients - usual scalar 
product-like formulas used for initial data and the modes 
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lead to divergencies Q, El • 

Second problem is completeness of the sum of quasinor- 
mal modes, i.e. if every solution of the wave equation can 
be (uniquely) represented by QNMs. It turns out that it 
essentially depends on the analytic and asymptotic prop- 
erties of the potential, e.g. long range potentials usually 
make QNMs incomplete (by causing a tail due to a scat- 
tering at large distances [l4L Il5jl and discontinuities in 
the potential improve completeness (by increase in the 
number of modes [Tri H^P- 

In this Article we address the first problem and not 
the second. In particular, we show that the divergen- 
cies have their origin in too careless use of analogies with 
normal modes, thus leading to mathematical inconsis- 
tencies. In the evolution of smooth initial data, in pres- 
ence of a well behaving potential, divergencies are not 
expected. We perform more careful calculations and ob- 
tain a well-defined representation of solutions, which in 
some region of spacetime (roughly, common future light 
cone of all initial data) takes the form of a sum over 
spatially truncated quasinormal modes with constant co- 
efficients. The coefficients can be explicitly calculated 
from the initial data by a scalar product-like formula in- 
volving initial data and QNMs. Though in general, they 
may depend on time and space, but it turns out that they 
become constant in a special region of spacetime, where 
the quasinormal modes become meaningful. The appear- 
ance of spatially truncated profiles of QNMs solves the 
problem of divergence at large distances and restores the 
principle of a causal propagation. We argue that this is 
the correct way to unterstand the quasinormal modes (cf. 
0) since the result is a natural consequence of solving 
the wave equation Q via Laplace transform and treating 
the asymptotic behaviour of the integrand in the inverse 
transform with appropriate care, with no other ad hoc 
assumptions about the form of the result. 

The problem of completeness 0,^3 wm * appear in the 
calculations. We will introduce necessary assumptions 
to eliminate it. As will be more transparent during the 
calculations, these two problems are quite unrelated and 
therefore we find it justified to study them separately. 

Considerations presented in this Article are not fully 
rigorous. We are still working on a corresponding theo- 
rem, which shall contain analysis of the class of evolved 
initial data (domains of operators) and conditions on 
the analytic and asymptotic properties of the potential. 
Therefore we would characterize these results as one step 
more rigorous than other general results known so far. 

This paper is organized as follows. In the second sec- 
tion we revise the theory of normal modes in order to un- 
derstand (later) some essentiall differences between nor- 
mal and quasinormal ones which might easily be over- 
looked in a less rigorous approach. In the third section we 
present a general theory of quasinormal modes and em- 
phasize restrictions in its validity. In the fourth section 
we present our main result - calculation of the Green's 
function, which contains construction of its asymptotic 
form for big s (Laplace transform parameter) and appear- 



ance of quasinormal modes with constant coefficients in 
some region of a spacetime. In the fifth section we briefly 
discuss a connection between the region of validity of the 
quasinormal mode expansion, its completeness and prop- 
erties of the potential. In the sixth section we present an 
example, which can be solved analytically, with results 
of a numerical simulation, which support our theoretical 
picture. 



II. NORMAL MODES 

The simplest proofs of normal mode expansions rely on 
the properties of the corresponding operator —d% + V(x) 
defined together with boundary conditions on a given 
space of functions. If the operator is self-adjoint the 
spectral decomposition into its eigenfunctions leads im- 
mediately to the normal mode expansion. Since in the 
case of quasinormal modes the corresponding operator 
is no more self-adjoined (different boundary conditions) 
no simple eigenfunction expansion exists and hence one 
cannot take much advantage of the 'algebraic' methods. 
This is the reason why we present an 'analytic' proof 
of normal mode expansion. The same technique will be 
later used for quasinormal modes. 



A. V = 

Consider the wave equation in one dimension on an 
interval [a, b) 

i){t,x)-il)"{t,x)=Q. (3) 

with Sturm-Liouville type boundary conditions at x = a 
and x — b and initial conditions: -0(0, x) = f(x), 
tp(0, x) = g{x). It is well known that the solution can 
be expressed as a sum over normal modes 

iP{t,x)=J2a n Mx)e luJnt , (4) 

n 

where <fi n (x) are normal modes, solving a corresponding 
eigenvalue problem (obtained from by time separa- 
tion) with an eigenvalue w„, which plays a role of a (real) 
frequency of the mode. a n are excitation coefficients, 
uniquely determined by the initial data f{x),g(x). Al- 
though the result is standard, we want to derive it here, 
in order to be able to compare it to the derivation of a 
quasinormal mode expansion and understand the differ- 
ences. That is why the method chosen here is the same 
as the one introduced by Nollert and Schmidt Q for def- 
inition of quasinormal modes. 

Starting from @ we perform a Laplace transform of 
the wave function 

/>oo 

4>(s,x) = £[tp(-,x)](s) = / e~ st ip(t,x) dt (5) 

Jo 
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and obtain an inhomogeneous ODE 

s 2 ip(s, x) - $"(s, x) = sf(x) + g(x) = j(s, x) (6) 

This equation can be solved by a Green's function tech- 
nique 

[s 2 -d*]G(s,x,x') = S{x-x'), (7) 



ip(s,x) = / G(s,x,x') j(s,x') dx' . (8) 



A proper Green's function has to be constructed, i.e. 
which, beside the differential equation 0, must also sat- 
isfy the same Sturm-Liouville boundary conditions as tj) 
does. As is well known, too, the following construction 
leads to the required result Take two solutions of 

a corresponding to © homogeneous ODE (with j(s,x) 
replaced by zero) 



[ s 2 -dl] Ms, 



o, 



which satisfy one boundary condition each, </>_ the left 
one (at x = a) and <j) + the right one (at x = b). Then 
the proper Green's function reads 



G(s, x, x') 



</>_(s,x<) <^ + (s,a: > ) 
W[i-J+](s) ' 



(10) 



with a;< = min(x,2:') and x > = max(x,a;'), and 
W[$-{x),$+(x)](a) = 0_(ar) $' + (x) - $'_(x) $+(x) a 
Wronski determinant, independent of x. 

Now, we are ready to invert the Laplace transform 



4>(t,x) =C- l ${;X)]{t) 



1 



e st tp(s,x) ds, (11) 



where the integration in the complex plane is to be per- 
formed parallel to the imaginary axis, to the right of all 
poles (see Fig. [IJ. Making use of the relation JSJ and the 
definition of j(s, x) we can express it, after some manip- 
ulations, as 



ip(t,x) 



d_ 

dt 



C- l [G{;x,x')]{t) f{x') dx 1 



(12) 



or shortly 



£- x [G{;X,x')}(t)g(x') dx', 



-L-^G]* f + L-^G}* g 
^-G * / + G * g. 



(13) 



Here, the initial data separate from the rest of the prob- 
lem, i.e. the boundary conditions and the wave equation 



lm s 



Re s 



FIG. 1: Contour of integration (C) for the inverse Laplace 
transform. 



itself. The whole universal information about solutions 
of the problem is contained in the Green's function 

1 



(9) G{t,x,x') = C- x [G{;x,x'))(t) 



2ni 



c 



e st G(s,x,x') ds 
(14) 



or, more explicitly, 

G(tl x, *') = —/ e°< *-(*,*<) *+{*,*>) ds _ (15) 
2m Jc W[<t>-,(j>+](s) 

A theorem of Weyl says that the functions <j>±(s,x) 
are entire functions of the parameter s. Therefore the 
function G(s,x, x') and the integrand are meromorphic 
in the whole complex plane, with possible poles only at 
zeros of the Wronskian W [(/>-, </>+](s). 

The simplest way to calculate the integral is to close 
the contour by a semicircle at infinity, either at the left 
or right hand side of the complex plane (Fig. [IJ. The 
contribution of the integral along the semicircle depends, 
in general, on the values of t,x and x' . Indeed, it can 
be seen by an example, where the boundary conditions 
ip(t, a) — ip(t, b) — have been chosenjS], that 



G(s, x, x) 



= sO<-a) 



s(fe-x>)l 



2 s [ e s(a-b) _ e s(6-a)] 



(16) 



(more details can be found in Appendix |A*J . Its asymp- 
totic form for |s| 3> 1, Re(s) < is 

p s\x— X \ 

G(s,x,x')~ 77 — (17) 



2s 



and for \s\ > 1, Re(s) > 
G(s, x, x) 



2s 



(18) 



It implies that the integrand in (|15fl behaves asymptoti- 
cally like 



a s(t+\x— x'\) 



2s 



-, Re(s) < and 



„s{t—\x—x'\) 

2s • 



Re{s) > 0, 
(19) 
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respectively. The first function, for Re(s) < 0, is always 
small for t > and the integral over the left semicircle 
vanishes for all t > 0. The second function, for Re(s) > 0, 
is small only for t < \x— x'\ and the integral over the right 
semicircle vanishes only for t < \x — x'\. 

Since, by definition of the inverse Laplace transform, 
all poles in l|15|) lay to the left of the original contour C, 
closing the contour to the right is free from the contribu- 
tion of poles. In contrary, closing the contour to the left 
produces an additional contribution of residua calculated 
at the positions of the poles. The residua correspond to 
the values of s at which the Wronskian W[<fi-, 4>+}{s) van- 
ishes. At these points both functions are proportional 



<j) + (s n ,x) = Cn<f)-(s n ,x) = C n (j)n{x) 



(20) 



with some proportionality constant c n . The integral 
in l|15fl over a closed (to the left) contour may, by the 
Cauchy theorem, be expressed as a sum over the residua 



Introducing the excitation coefficients 



A n = a n \ [s n f(x') + g(x')] <j) n {x) dx' (26) 

J a 

the solution takes a compact form 

iP(t,x) = J2 A n<f>n(x)e Snt - (27) 



This result shows that (i) the form of the modes de- 
pends only on the form of the wave equation and bound- 
ary conditions and not on the initial data, (ii) the excita- 
tion coefficients are uniquely determined from the initial 
data, (iii) the expansion in a sum over modes is possi- 
ble for every time t > 0, and (iv) the causality condition 
holds, i.e. initial data may influence the solution only 
within their future light cone (cf. 1241) ). 



G(t,x,x') = 

— • 2iri Res 



2m 



) <j>-(s,x<) ^+(s,a: > ) 
W(s) 



^e s " t ^_(s„,x < ) (j> + (s n , X>)Res 

n 

^V"*c„ <£„(£<) <t>n( x >) w ^ s y 



W{s) 



Introducing "normalization constants" a Tl 
— c n /W' {s n ) we finally arrive at 

G(t,x,x') = *^2a n e Snt (/) n (x) (f> n {x') 



(21) 



(22) 



for every t > 0, while closing the contour to the right 
gives simple 



G{t,x,x') = 



(23) 



for t < \x — x'\. This is the main result for normal mode 
expansion. It is important to note, that for t < \x — x'\ 
the contour me be closed to either side, giving either a 
sum over modes or zero in the result, what leads to the 
conclusion that the sum over the modes converges in this 
case to zero. Combining both results into one formula 
gives 

G(t,x,x') = eit-lx-x'D-J^ane^Mx) Mx'), (24) 

n 

making the causality property explicit. 

For completeness, let's calculate, in the last step, the 
wave function, using (|12|) 



/ ^a n s n e Snt (j> n (x) (f) n (x') f(x') dx' 
■> a n 
,6 

+ / ^2a n e Snt (f> n {x) (f>n(x') g(x') dx' . 



(25) 



B. V + 

An analogous construction can be carried out if the 
waves propagate in presence of a potential, like in Q, 
as long as the boundary conditions are of the Sturm- 
Liouville type and are imposed at the ends of a finite 
interval. 

Let the boundary conditions have the general form: 



/3 ip(t,a) - a ^-(t,a) = 
ox 

S^(t,b)-j^(t,b) = 0. 



(28) 
(29) 



with a,P, j, 5 G R. By the technique described in more 
detail in section IV it can be shown that 



4>-{s,x) 



-s(x—a) 



{l + Ois- 1 )) 



4>+{s,x) 



+ e s{x-b) I 



(l + Ois- 1 )). 

(30) 



Then, first the approximate Wronskian W[(/)-,<j>+](s) 
and next the approximate Green's function can be ob- 
tained with the same precision (see Appendix El f° r 
more details). Considering its asymptotic behaviour for 
\s\ 1 one finds the same formulas as in the previous sec- 
tion (V = 0), namely (UJ and JTHJl. Their consequence 
is the normal mode decomposition 



G(t,x,x') 



4>n{x') 



(31) 
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for every t > and 

G{t,x,x')=0 (32) 

for t < \x — x'\. For nonnegative potentials all s n must be 
purely imaginary (because — s„, the eigenvalues of —d 2 + 
V(x), must be real and positive). For negative potentials 
s n must be either purely imaginary or real positive (in 
all cases s 2 n < — inf xe [ ajb ] V(x)). 



for every t > and 

G(t,x,x')=0 (39) 

for t < \x — x'\. The essential difference is that the fre- 
quencies s n are no longer purely imaginary or real. They 
can be complex and for nonnegative potentials must lay 
in the left half of the complex plane (Re(s n ) < 0). In case 
of positive potentials there may be additional s n which 
lay on the positive half of the real axis (s ra € R + ). 



C. Leaky cavity 

Even for some non-Sturm-Liouville boundary condi- 
tions, which lead to non-self-adjointness of —d 2 + V(x) 
on the corresponding domain, the mode expansion may 
be possible for all times t > 0. Again the same tech- 
nique of constructing approximate Green's function can 
be used for the proof. 

Consider boundary conditions of the form 



iK*.o) = o, 



dip 
~dt 



(33) 
(34) 



The left boundary x = is fully reflecting and the right 
one is leaky - part of the incoming signal from the left 
will be reflected and part will be transmitted and lost 
from the domain [a, b]. The parameter 7 measures the 
leakage: for 7 = the boundary is fully reflecting and 
for 7 = 1 it is transparent for waves coming from the left. 

The boundary conditions for ip(t, s) translate into con- 
ditions on the Laplace transformed functions 

0_(a,O) = 0, (35) 

s4> + ( s ,L)+ 7 ^±(s,L) = (36) 

(we have assumed ip(0, L) = for simplicity, otherwise 
this term should appear on the r.h.s. of l|C3J) ). 
The solutions are 

4>-(s,x) =e sx -e~ sx (37) 



<i>+(s,x) 



f{x-L) ( 7 _ 1) + e -»(«--t) ( 7 + 1) 



Next, the Wronskian W[(/>_,</>+](s) and finally the 
Green's function can be constructed (see Appendix IClfor 
more details). The asymptotic behaviour for |s| 3> 1 is 
the same as for normal modes discussed above and given 
by (|T7jl and ijTBl) . Analogously the mode decomposition 
follows 



G(t,x,x') 



a n e Snt cj> n (x) 4> n {x') 



(38) 



III. QUASINORMAL MODES 

Consider now the same wave equation with a potential 



d 2 



if 



dt 2 dx 2 



ip(t,x)=0, 



(40) 



but on an infinite axis, and initial conditions 

M0,x) = f(x), i>{0,x)=g(x). 



(41) 



About the potential we assume that it is real, continuous 
and vanishes faster than l/\x\ as x — > ±00. To construct 
a similar set of modes a kind of a boundary condition at 
x —> 00 is needed. Mathematically it corresponds to the 
definition of a domain of solutions. Historically, one has 
used a boundary condition of an asymptotically outgoing 
wave. This bases on the experience from the scattering 
theory. Although it is to some extent correct, it proves 
unsatisfactory in more detailed calculations, namely, for 
Re(s) < this boundary condition fails to be unique. 
Nollert and Schmidt proposed a procedure, based on 
a Laplace transform and analytic properties of the trans- 
formed Green's function, which leads to uniquely defined 
quasinormal modes. We have demonstrated the method 
of a Laplace transform in the previous section. What 
remains, is to see how the change of the boundary condi- 
tions (and therefore the domain of solutions) causes the 
quasinormal modes to replace the normal ones. 

Performing again the Laplace transform, which is de- 
fined for Re(s) > 0, 



$(s,x) = C[ip(-,x)](s) 



we obtain the inhomogeneous ODE 



ip(t, x) dt, (42) 



[-81 + s 2 + V(x)] ^(s, x) = sf(x) + g(x) = j(s, x). 

(43) 

We solve it by a Green's function technique 

+ s 2 + V(x)] G(s, x, x') = 6{x ~ x'), (44) 



ip(s,x) = / G{s,x,x') j(s,x') dx' . 



(45) 
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The proper Green's function, again, is constructed from 
two solutions of the corresponding homogeneous ODE 



[-dl + s 2 + V(x)] 4>±{s,x) =0, 



(46) 



which have the property, that 4>±(s,x) satisfies the 
boundary condition as x — > ±oo, respectively. The 
boundary condition is chosen to be the square integra- 
bility at infinity, i.e. 



\4>±(s, x)\ 2 dx < oo 



(47) 



for every a £ R. As long as Re(s) > 0, it is equivalent to 
the outgoing wave condition. It can be seen by an explicit 
construction of the asymptotic behaviour of the solution 
which is square integrable at the given end. Since the 
potential vanishes at infinity, we find 



4>-(s,x) 
4>+(s,x) ~ e- sx 



as 
as 



x 

X 



-oo 
-oo. 



(48) 
(49) 



The inverse Laplace transform restores the time depen- 
dence, roughly, in the form of a multiplication by exp(si) 
(in the case of a single mode, at least), what gives 



(j>-{t,x) ~ e s{t+x) 
<Mi,x)~e^-*) 



— oo 

+O0. 



(50) 
(51) 



Both are purely outgoing waves, what agrees with the 
original idea of the quasinormal modes. 

A comment on a difference appearing in (|46|l relative to 
the previous section is in place here. Although the equa- 
tion is nearly identical with the differential operator 
—d 2 + V(x) is no more self-adjoint on the domain defined 
by the outgoing wave boundary conditions. Therefore its 
eigenvalues —s 2 do not have to be real. For normal modes 
it holds s — iuj, to G M, what makes them purely oscilla- 
tory in time with the frequency lo. This behaviour is not 
expected by quasinormal modes. 

By these boundary conditions, the Green's function 



G(s, x, x') = — 



<j>-{s,x<) 0+(s,x>) 
W[4-J+](s) 



(52) 



with a;< = min(a;, x') and x > = max(i, x'), is square in- 
tegrable in both variables x and x' . So is the transformed 
wave function ip(s,x) calculated from Ij45(l if j(s,x) (or 
equivalently f{x) and g(x)) is of compact support. 

The inversion of the Laplace transform goes the same 
way as in the preceding section, following exactly the 
formulas (| 1 1 11 - (| 1 4|) . and leads to the integral 



G(t,x,x') = 



-1 



■ (S,X<) <Ms,3b 

W[$-,$+](s) 



ds. (53) 



Here, the same idea of closing the integration contour by 
a semicircle at infinity seems most natural, but in this 



case the dependence of (f>±(s,x) on s is much less trivial 
than it was by normal modes. First, there is no theorem 
saying that they are entire functions of s. The special 
boundary conditions make them more complicated func- 
tions and poles as well as essential singularities producing 
branch cuts may occur. Fortunately, the poles in </>±(s, x) 
cancel with those in the Wronskian W[4>-, cf)+](s), but es- 
sential singularities and cuts in general survive. 

Second, it is not easy to determine the asymptotics of 
the integrand in (|53J) at big \s\. Its control is essential 
for estimation of the contribution of the integral along 
the semicircle at infinity. Let's first shortly recall what is 
popularly written in the literature: reviews the state 
of the art in the theory of quasinormal modes. It con- 
tains discussion about the divergence problem and spec- 
ulations about time-dependence of the excitation coeffi- 
cients (proposed by N. Andersson [13] )• Though, it does 
not bring firm answers, since on the mathematical level 
the problem is not treated carefully enough. 

It looks natural that at positive times t > the ex- 
ponent exp(si) in H53J) goes to zero for big \s\ with 
Re(s) < 0, so if the rest of the integrand is somehow 
bounded (here is the point!) then the contour can be 
closed by a semicircle at infinity on the left hand side. 
An important point is that the integrand is defined, via 
the Laplace transform, only for Re(s) > and its do- 
main of definition must be extended to the left half- 
plane Re(s) < 0. It can be done in various ways, al- 
though e.g. imposing the condition of square integrability 
proves wrong and useless, and the one of the asymptotics 
cj)± ~ exp(=psa;) at x — > ±oo seems to be mathematically 
not well-defined, because the other linearly independent 
asymptotic solution exp(±sx) is subdominant relative to 
it and cannot be eliminated |25j . The only correct method 
is the analytic continuation of 4>±(s, x) from the right to 
the left half of the complex plane of s. Only this allows to 
use the Cauchy theorem and replace the complex integral 
by a sum over residua. 

Now, assume for a moment that there are no cuts in the 
complex plane, and use the Cauchy theorem. We arrive 
at the sum over modes, like in the previous section, 



G(t,x,x') = 



but now the functions <j) n {x) are quasinormal modes. 
They are also defined by 4>±(s, x) calculated at the posi- 
tion of the poles s — s n in the left complex half- plane, but 
they do not satisfy the boundary conditions in their orig- 
inal form of square integrability. Since we have chosen 
analytic continuation, their asymptotic form as x — > ±oo 
for Re(s) < remains explicitly the same (analytically 
continued, too) as it is for Re(s) > 0. Moreover, every 
mode satisfies both boundary conditions, what gives 



t (x — > -oo) 



and <j) n {x — * +oo) ~ e 



(55) 

Since Re(s n ) < 0, the quasinormal modes are unbounded 
at both ends. Optically, it looks like the asymptotic 
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boundary conditions, criticized above as non-unique, due 
to an uncontrollable admixture of the other, subdomi- 
nant asymptotic solution. Here, however, thanks to the 
analytic continuation from the right half-plane, where 
the asymptotic conditions are unique, the construction is 
unique, too. The spatial unboundedness throws a shadow 
of a doubt on their applicability in the expansion of so- 
lutions of the wave equation. Before we come to that 
question, let's first observe that already the assumption 
of boundedness of the "rest" of the integrand in H53[) is 
wrong. What we can expect, is rather a boundedness by 
an exponent ~ exp(±sa:). Therefore, the procedure of 
closing the contour on the left and assuming that the con- 
tribution of the semicircle integral vanishes, seems only 
reasonable for late times t > \x — x'\. 

To answer the question more exactly, for which values 
of t, x, x' the contour can be closed left or right with 
finite or vanishing contribution from the corresponding 
semicircle integral, we have to analyze the asymptotics 
of the integrand at big \s\ in more details. 



IV. CONSTRUCTION OF THE GREEN'S 
FUNCTION 

An essential observation is that the question if the con- 
tour of integration in the inverse Laplace transform of the 
Green's function may be closed on the right or left de- 
pends only on the asymptotics of the Green's function 
for big \s\. This, further, depends on the asymptotic be- 
haviour of <p±(s, x) for big |s|, which satisfy the ordinary 
differential equation 



d 2 J{s,x) = [s 2 + V(x)] $(s,x). 



(56) 



By means of the Liouville-Green approximation |21| we 
are able to construct the asymptotic solutions of this 
equation to the leading order in \s\ and control the er- 
ror. The same will be possible, in the next step, for the 
Green's function. 

The theory says that there are two, so called, dominant 
solutions 

Ms,*) = u- 1 / i -e+f ul/2dx ■(l + eifoa;)) (57) 
y 2 {s,x) = u^ A -e-^ 1/2dx -(l + e 2 ( S ,x)) (58) 

with the bounds on the error functions e\ , e 2 and their 
derivatives 



|ei(*,aO|, lu-^flUifoaOl < e T <^ - 1 
\e 2 (s,x)\, |ii- 1/2 a x e 2 ( S ,a;)| < e r "-' - 1 

with T a j) defined (for our purpose) by 

rb 



\F'(x)\ dx 



(59) 
(60) 



(61) 



and 



F(x) = / \u- 1/4 (u- 1/4 )" - vu- 1/2 ] dx. (62) 



The auxiliary continuous functions u and v can be chosen 
freely and in our case must satisfy the condition u(x) + 
v(x) — s 2 + V(x). The boundaries are a\ = — oo and a 2 — 
oo. The function Re J v}l 2 dx must be nondecreasing on 
the intervals (a%,x) and (x,a 2 ). T aitX must be finite for 
every x (i = 1, 2). 

In our problem the most convenient choice is: 



u(x) — s 2 , v{x) — V(x), 



(63) 



with the choice of the complex root it 1 / 2 = s. Then the 
solutions take the form 

<Pi(s,x) = 4= e +ax -(l+ei (*,*)) (64) 

V2(s,x) = -±=e- sx -(l + e 2 (s,x)), (65) 

The function Re J u 1 ^ 2 dx = Re(s)x is obviously nonde- 
creasing as x increases on (— oo, oo) as long as Re(s) > 0. 
The auxiliary functions for the error estimates take quite 
simple form, too, 



F\x) 



V(x) 



(66) 



T 0l , x = T^oo,, = t~t f \V(x)\ dx (67) 

\S\ J-oo 
1 f + °° 

Ta 2 ,x = T +00 . x = — / \V(x)\ dx, (68) 
\s\ J x 

and the errors are bounded according to f59"|l . JHOI)- Here 
it becomes transparent, that this construction works only 
for integrable potentials V £ L 1 (R). 

It can be easily seen from the bounds on £1,62, that 
the solutions ipi , <p 2 behave asymptotically like 



(fi(s,x) 
ip 2 (s,x) ~ 



as x — > —00 
as x —i- +00 



(69) 
(70) 



and therefore satisfy the boundary conditions. Hence, we 
choose 

4>-(s,x) = <pi(s,x) 4>+(s,x) = (p 2 (s,x). (71) 

Now, the behaviour of the solutions 4>± at big \s\ may be 
found. Since the potential V(x) is continuous and van- 
ishes asymptotically like | V(ac) j < C with some 
C,e > for \x\ -> 00, the integral ||V|| = f+™ \V{x)\ dx 
is finite. We obtain the following estimates 



T < 



T < 

^ +00, X — 



(72) 



and the bounds on the errors 



| £l ( s ,x)|, lu-Va^^a;)! < ell v H/l-l-l (73) 
|e 2 (a,a;)|, \u-^ 2 d x e 2 (s,x)\ < e ll^ll/l-l _ 1. (74) 
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These can be expanded for \s\ 3> 1 to give 
|ei(s,x)|, lu-^d^x)] = o(^ 
\e 2 (s,x)\, \u-^ 2 d x e 2 ( S ,x)\ = O 



Therefore 

4>-(s,x) 

4>+(s,x) 



1 
1 

7^ 



1 + 



1 + 



(75) 
(76) 

(77) 
(78) 



Further, it can be shown, that the error of their Wron- 
skian is bounded in a similar way 



Wl 



,4>+](8) = -2- 



l + O 



(79) 



Finally, the asymptotic form of the Green's function for 
big \s\ may be written 



G(s, x, x') — 



-s\x— x'\ 
~2~S 



l + O 



(80) 



It is important to note, that it asymptotically, for |s| — > 
oo, coincides with the free Green's function, i.e. for a 
wave equation with V{x) = 0. It is not unexpected, 
since in the r.h.s of Ij56(l . in s 2 + V(x), the potential V(x) 
is small compared to s 2 as \s\ — > oo (more exactly s 2 3> 



Unfortunately, this approximation is valid only for 
Re(s) > 0, as already noted. The choice of Re(s) < in 
the construction of the dominant solutions leads to an ex- 
change of the formulas for ipi, (p 2 . The reason is 
that the function Re J u 1 / 2 dx must be nondecreasing in 
x, so it requires a change of u 1 / 2 from s to —s. The whole 
construction can be repeated with the exchanged formu- 
las, but then the functions <j)± are constructed separately 
in the left and right complex half-plane of s. Obviously, 
they are not analytic in s, or more precisely, the func- 
tions (f>± constructed for Re(s) < are not the analytic 
continuations of those constructed for Re(s) > 0, what 
is required. 

Having this in mind, let's come back to the analysis 
of the inverse Laplace transform for the Green's function 
The integrand equals approximately, for big \s\, 



„s{t—\x — x'\) 

2s 



l + O 



(81) 



what means, that the contour of integration may be 
closed on the right for t — \x — x'\ < 0. That leads to 
the conclusion 

G(t, x, x') = for Q<t<\x-x'\. (82) 

It is nothing else than the causality condition. The same 
result is obtained immediately in the free case (V(x) = 




FIG. 2: Validity region (shaded) for the quasinormal mode 
expansion. 



0) . It is not surprising that the presence of a potential in 
the wave equation does not change the causality of the 
wave propagation. 

Until now, we have reproduced "one half" of the re- 
sult holding for normal modes (|23|l . namely the causality 
condition for the Green's function. The "second half" 
(eq. (|22|l for normal modes) should be the expansion of 
the Green's function into quasinormal modes. Unfortu- 
nately, the situation here is much more difficult: not only 
an expansion for any t > cannot be shown, because the 
integral along the semicircle at infinity on the left hand 
side may diverge, but also the expansion for t > \x — x'\, 
where the Green's function is nonzero, cannot be shown 
always. It can be proved only under further assumptions 
on the potential V . In general, it is expected that the 
expansion holds for some t > t (x,x') > \x — x'\. In 
the Poeschl- Teller potential has been studied in detail 
and the function to(x,x'), which gives the earliest time 
of validity of the expansion, has been found. This region 
(i, x) in evolution of initial data with compact support 
has the form like shown on the Fig. [21 

We want briefly sketch the technical difficulties which 
arise. The functions <p±{s,x), which for a given s solve 
the ordinary differential equation l|56|l with one bound- 
ary condition, at x — > ±oo, respectively, have a simple 
asymptotic behaviour on that end 



4>-(s,x) 
(f>+{s,x) 



+oo, 



but a complicated behaviour on the other end 



<p-(s,x) - A{s)e sx + B{s)e- S 
j>+(s,x) - C(s)e- sx + D(s)e s 



x ■ 

X 



(83) 
(84) 



-oo (85) 
-oo. (86) 



The "transmission" coefficients A(s) and C(s) tend to 
1 and the "reflection" coefficients B(s) and D(s) tend 
to zero as \s\ — > oo, but the rate of their convergence 
(~ s~ n ) is usually too slow to make the corresponding 



terms appearing in the Green's function, like 

B{s)C(s)e- s{x+x,) or A(s)D(s)e s( - x+x '\ (87) 

to vanish as \s\ — > oo. This corresponds to the scattering 
of the waves on the potential, even in the limit of high 
frequencies. It causes that the contour of integration can- 
not be closed on the left hand side immediately after t 
passes | a; — x'\, what is the case for a free wave equation 
(exp[s(t— \x — x'\)] — > 0), but only later, when other con- 
tributions to the Green's function (~ exp[s(t — x — x')] 
or exp[s(t + x + x')]) vanish, too. Finally, the contour 
may be closed, and by the same, the expansion is valid, 
for times t > to(x,x') > \x — x'\ (Fig. [2J. The function 
to(x,x') is in general a complicated function of the po- 
sitions x, x' and depends essentially on the shape of the 
potential V(x). 

Another conclusion can be drawn from this construc- 
tion. If the contour can be closed on the left hand side 
and the expansion into quasinormal modes holds, the 
modes will only appear in a spatially truncated form. 
The reason for the truncation is the causality condition 
together with the assumption that the initial data have 
compact support. They can influence only the region of 
spacetime laying inside their future light cone. In the 
outer region the solution must remain zero. In an op- 
timal case (which can be realized by some potentials), 
when the expansion holds for all t > to(x,x') = \x — x'\, 
the Green's function can be written as 

G(t,x,x') = 6(t-\x-x'\)-J2<Xne Snt <f> n (x) M*')- ( 88 ) 

n 

It is the same result as for the normal mode expansion, 
but here the 9— function plays an additional role: it trun- 
cates the modes at finite times and avoids their diver- 
gence at big distances (Fig. |3J). Furthermore, it is not 
indifferent now if we write the 9— function in front of the 
sum or not, as it was by normal modes, because the sum 
itself is divergent just there, where 9{t — \x — x'\) equals 
zero, i.e. for t < \x — x'\. The sum over normal modes, 
in contrary, converged to zero in that region. 

The solution of the wave equation can be written in 
the form 

1>(t, x)=J2 B n{t, as) Mx) e s "*, (89) 

71 

which is similar to that for normal modes, but the ex- 
citation coefficients become in general time and space 
dependent and are defined by 

pX + t 

B n (t,x) = a n / [s n f(x') +g{x')\ 4> n (x') dx' 

Jx-t 

+ a n [f(x - t) cj) n {x - t) + f(x + t) cj) n {x + t)] , 

(90) 

what results from the presence of the 9— function in (18811 . 
The B n (t,x) coefficients are zero outside the future light 
cone of the initial data. They get constant at points in 




FIG. 3: Appearance of truncated quasinormal modes in evo- 
lution of initial data with compact support. 



spacetime (t, x), where all initial data lay inside the past 
light cone of (t, x), i.e. when all initial data "can be seen" 
(Fig. E| . 



A. Case V = 

It should be noted that the general result l|88|) contains 
the case of a free wave equation with V(x) = 0. In this 
situation the exact solutions of (|56|l are 

h(s,x) = e* sx (91) 
W$-,$+]{a) = -2s (92) 

(where we have multiplied both solutions by y/s for con- 
venience, what changes the Wronskian accordingly but 
has no influence on the Green's function). Then the 
Laplace transformed Green's function equals exactly 

G(s,x,x')= e SW Xl . (93) 
2s 

The only pole is placed at s = 0. The corresponding 
quasinormal mode is trivial 

Mx)=l, (94) 

but what makes the solution of the wave equation non- 
trivial is the #(...) factor in front of the QNM sum 

G{t,x,x') = 9(t-\x-x'\)-a Q e Sot <j)o(x)Mx') 

= \ 6{t-\x-x% (95) 

where «o = Res[l/VF(s)] s=Sf) = 1/2. It is the well known 
Green's function for a free wave equation, which cannot 
be reproduced by formula l|54|) valid only for long times 
t > \x-x'\. 
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V. VALIDITY OF THE EXPANSION 

A discussion when the quasinormal mode decomposi- 
tion is a full solution of the Cauchy problem is in place 
here. Beside the requirement of compact support of 
the initial data, in question are conditions on the po- 
tential V(x). In this construction it has been first as- 
sumed that V must be continuous and integrable, what 
is not a strong restriction. Though, further it has been 
assumed that the analytic continuation of the Laplace 
transformed Greens's function G(s, x, x') should have no 
complex branch cuts. It is a complicated indirect con- 
dition on the potential. We do not want to go into a 
detailed discussion of this problem in this short Article. 
Its role is rather to clarify some questions regarding QNM 
decomposition on a heuristic level. We are working on a 
result in form of a theorem where the problem of branch 
cuts will be addressed with adequate care. 

However, some general conclusions can be drawn with- 
out much analytical work. It has been observed that 
smooth potentials with exponential fall-off at infinity pro- 
duce no cuts in G(s). On the contrary, potentials with 
power-law decay usually produce a complex cut emanat- 
ing from an essential singularity at s = 0. This relation 
can be explained by the result of Hod [l4| which relates 
the fall-off of the solution at late times with the asymp- 
totic behaviour of the potential (see also 0). Roughly, 
initial data evolving in presence of exponential poten- 
tials decay exponentially (like in Poschl- Teller potential 
studied in [t|) and in power-law potentials decay with 
a power-law (what is well known for black-hole back- 
grounds - see |a] for detailed analysis of the Schwarzschild 
case and [2^ for a review of all other black-holes). Since 
a sum of QNMs decays exponentially in time while the 
practise shows that solutions with a branch cut contri- 
bution in G(s) fall-off slower, like power-law, it becomes 
plausible that exponential potentials should not give rise 
to a cut while long-range potentials produce brach cuts 
necessarily. 

The second open problem is that of the earliest time 
to(x, x 1 ) > \x — x'\ when the QNM expansion begins to be 
valid. Its dependence on the shape of the potential seems 
to be very complicated, depending probably strongly on 
the analytic properties (differentiability) of V(x). It is 
related to the scattering of waves in the limit of very 
high frequencies, where hopefully some results from the 
scattering theory may be used. It is known that discon- 
tinuities in the potential or its derivatives increase the 
number of modes 0, 0] . 



VI. EXAMPLE: V = Voe~ 2r 

In this section we present an example of a wave equa- 
tion, which can be, to a large extent, solved analytically. 
We choose 



and modify a little the original setup of the problem. 
Let now the spatial coordinate be defined on r € (0, oo) 
and let us choose the boundary condition at r = to be 
ip(t,0) = 0. Such conditions correspond, for instance, to 
a 3-dimensional wave equation with a spherically sym- 
metric potential (after separation of the angular depen- 
dence and multiplication of the wave function by r). In 
the 1-dimensional problem on the half-axis the boundary 
at r = acts reflecting on the incoming waves. 

This potential has the nice property that it does not 
produce any essential singularity, and thus no cuts, in 
the Green's function. The reason is its fast decay at 
infinity, i.e. at the open end. Therefore the solution 
can be expanded purely into quasinormal modes, without 
contribution of an integral along the branch cut, what 
allows a better insight into the problem of completeness 
of the expansion. Another reason is that the asymptotic 
behaviour of solutions <fi ± (s, x) for big s is relatively 
simple. 



A. Analytical results 

First observation is that only one end, at r — > oo, is 
open in this example, while the other one at r = is 
reflecting. It leads to a modified free Green's function 

G (s, r, r') = 1 e''^ - 1 e"^'). (97) 
2s 2s 

The first term is the usual one for a free wave equation 
and the second one describes the reflection of waves at 
the origin. 

Following the way of constructing solutions described 
in the previous sections, we present here only results of 
the consecutive steps. The solutions of the homogeneous 
ODE, satisfying one of the boundary conditions each, can 
be expressed by means of the Bessel functions 



M*,r) = J M N s {qe' r ) - N s (q) J s {qe~ r ) 
h{s,r) = J s (qe- r ), (99) 

where Vq = —q 2 , i.e. the potential is negative p^. Their 
Wronskian equals 



W[<j> ,<P + ]{s) = -J s {q) 

7T 

and the Green's function is 

tt J_ s (qe- r <) J s {qe- r >) 



(100) 



G(s,r,r') 



2 sin(7rs) 
vr J- S (q) 



2 sin(7rs) J s (q) 



J s (qe r ) J s (qe r ). 



(101) 



It can be shown that 



V{r) = V e 



-2r 



(96) 



G(s,r,r') - 1 e- s l"'l - 1 e-< r +^ = G (s,r : r') 
2s 2s 

(102) 
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as \s\ — > oo. Quasinormal mode frequencies correspond 
to poles of the Green's function and these are at the 
values of s where J s (q) = 0. Therefore the quasinormal 
modes take a simple form 

Mr) =J a Jqe- r ). (103) 

They satisfy both boundary conditions 

</> rl (0) = and <j) n (r) ~ e~ SnT as r -> oo. (104) 

In this example, there exist additional poles at the points 
where sin(7rs) = 0. The residua calculated at these points 
cancel in integration of G(s, r, r'), but they do not cancel 
in integration of e st G{s,r,r '), which is to be performed 
for inverting the Laplace transform. They produce an 
additional series of modes, which shows up only at inter- 
mediate times \r — r'\ < t < r+r' (at earlier or later times 
they cancel) and contributes to the Green's function by 

\ ]T J n (qe- r )J n (qe- r ')e nt . (105) 

n— — oo 

This series converges in the given range of parameters 
t, r, r' . At r, r' 3> 1 it can be approximated and summed 
up to give a closed formula 

-~ + \h ( 9 e-(*+^')/ 2 ) + ij (qe«-r-r 'V 2 ) . (106) 

Finally, the Green's function is 

G(t,r,r') = 

'0, t<\r-r'\ 
1 

( - J] J n (qe~ r ) J n (qe~ r ') e nt , \r - r'\ < t < r + r' 

n— — oo 

J2a n J Sn {qe- r ) J Sn (qe~ r ') e s "*, i > r + r', 

(107) 

where a n = %N Sn {q)/ j^Js(q)\ s _ s ■ The solution takes 
the form 

4,{t,r) =Y / B n (t,r)J Sri (qe- r )e s " t 
n 

+oo (108) 
+ ^ C n (t,r)J n (qe- r )e nt 

TL— — OO 

with the excitation coefficients 



B n (t,r) = 




J Sn (qe- r ')[s n f(r')+g(r')]dr' (109) 



C n (t,r) = \fj^ Jn(qe- r ')[nf(r')+g(r')} dr' , (110) 

where we extended the definition of the initial data by 
f(r) = g(r) = for r < 0. 




FIG. 4: Appearance of truncated quasinormal modes in spher- 
ical symmetry. 



The formulas simplify if we assume we know the sup- 
port of the initial data explicitly 

supp(f) U supp(g) C [a, b}. (Ill) 

Then 

1. Zero: ip(t,r) = in two regions, t < a — r and 
t < r — b, which are causally separated from initial 
data. 

2. Sum over the additional modes: ^p(t, r) = 
^2,C n {t,r)J n (qe~ r )e nt for t < r + a, i.e. in the 
region, where no reflection at r = had occurred 
yet. 

C n become constant (t and r independent) in the 
region described by two conditions: r— a < t < r+a 
and t > b — r, i.e. where all initial data "can be 
seen" (contained in the past light cone) and before 
first reflected signal arrives. 

3. Sum over the quasinormal modes: ip(t, r) = 
B n (t,r)J Sn (qe~ r )e Srit in the region t > r + a, i.e. 
when some signal is already reflected at r = 0. 

B n become constant (t and r independent) for 
t > r + b, i.e. after all initial data have reflected at 
r = 0. 

The reflection at the boundary r = changes the pic- 
ture qualitatively (Fig. |2J . The quasinormal modes with 
constant excitation coefficients may first appear when all 
initial data reach the boundary and reflect. Therefore, 
the universal picture (sum over QNMs) shows up start- 
ing at the origin at time t — b and expands with the 
speed of light (here c = 1), creating a triangle (half light 
cone). After the time t, the QNM expansion is valid on 
re{0,t-b]. 
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B. Numerical simulation 

We have simulated this wave equation numerically. 
The initial data were chosen as 

{A ■ [1 - cos (fa;)] , for < r < 2, 

2 • A, ' for 2 < r < 22, 

A ■ [cos (f (x - 22)) + 1] , for 22 < r < 24, 
0, ~ for 24 < r, 

(112) 

with A = 1 and g(r) = 0. /(r) presents a long bump, 
starting from r — and extending to r = 24, smoothly 
connected with zero at both ends. On an interval [2, 22] 
its height is constant. The wave equation has been dis- 
cretized with grid steps dx = 0.02 and dt = 0.01, what 
guarantees stability of the numerical evolution. The spa- 
tial size of the grid was x £ [0, 100] (5000 grid points), 
what allowed a safe evolution until t = 76, when the first 
data reached the outer (free) boundary of the grid. 

The strength of the potential was chosen to be q = 2, 
what corresponds to Vq = —4. In this case (and in gen- 
eral, if Vq < 0) all poles are purely real. Numbering from 
right to left, i.e. from the dominant one (least damped) 
to more damped, their values are 

si = -0.2538, s 2 = -1.789, s 3 = -2.961, s 4 = -3.996, ... 

(113) 

It can be shown analytically, that the poles of the Green's 
function with big s (i.e. zeros of J s (q)) coincide approx- 
imately with the poles of T(s + 1), what gives s n — —n. 
Here, it agrees from n — 4 on within 1%. 

We have compared the results of numerical simulation 
with the solution I jlOSjl . We have calculated the exci- 
tation coefficients B n and C„ directly from initial data, 
taking into account only one, two or three least damped 
modes. It turned out that the differences in the damping 
exponents were so big, that the next modes, except the 
first one, were practically invisible. The reason was that 
when the QNM expansion began to be valid, it was t = 24 
and every mode was multiplied by exp(s n f). These fac- 
tors varied by ca. exp(— 24) from one mode to another, 
making practically impossible to observe any mode, ex- 
cept the first one. 

However, it is not to be seen as a problem. We want 
to stress that it is rather usual that only a small number 
of modes contributes essentially to the sum, while all 
other are damped fast to zero. It is a manifestation of 
universality and therefore it is a highly welcome property. 
Not to be groundless, we studied other values of q, like 
q = 2.4, for which si = -0.003126 and s 2 = -1.657. 
The first mode was so weakly damped, that it appeared 
as nearly static. Then, the dynamics was well described 
by the sum of the first and the second mode. 

Coming back to the case q — 2. Fig. [S] shows four 
"screen-shots" - shapes of the solution at a sequence of 
times (h = 26.71, t 2 = 29.77, t 3 = 33.52, t 4 = 38.86). 
The vertical segments show the right end of the light 
cone, i.e. r — t—b with b = 24, where the QNM expansion 
with constant coefficients is valid. Since the light cone 
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expands in time, the segment moves to the right. In the 
region r G [0, t n — b] we can see the solution approximated 
essentially by the first mode. Further to the right, a 
direct signal from the initial data, reflected at the origin, 
can be observed. Since it is expanded only formally into 
QNMs, with time and space dependent coefficients, its 
form is irregular and bears no sign of universality. 

Since the solution, as well as the first mode, are 
damped in time, they appear on every "shot" with a 
different amplitude. Dividing them all by the damping 
factor of the first mode, exp(sot n ), we can bring them 
to one scale, what is shown on Fig. |SJ The agreement 
is perfect and demonstrates that the solution is well ap- 
proximated by the first mode 

4>(t,r) SB, J Sl (qO e Sl * (114) 

with 

B\ = ai si / JsM^ r ') f(r') dr'. (115) 
Jo 
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VII. CONCLUSIONS 

The main result of this paper is the time-dependent 
Green's function l|88[l which allows, under further restric- 
tions on the potential V(x), to express the evolution of 
any complactly supported Cauchy initial data by means 
of a sum over spatially truncated quasinormal modes. 
The solution for the wave function takes the form l|89|) - 
(|9U[1 and is schematically shown on the figure [3J This 
clarifies the appearance of QNMs in the evolution at any 
time, solving the problem of divergence at big distances 
and restoring the principle of a causal propagation. The 
price is to allow the excitation coefficients depend on time 
and space in general, but it turns out that they become 
constant in a special region of spacetime (roughly, com- 
mon future light cone of all initial data) , where the quasi- 
normal modes become meaningful. We argue that this is 
the correct way to unterstand the quasinormal modes (cf. 
0) since the result is a natural consequence of solving 
the wave equation via Laplace transform and treating 
the asymptotic behaviour of the integrand in the inverse 
transform with appropriate care, with no other ad hoc 
assumptions about the form of the result. 

Unfortunately, the whole scheme is valid for all times 
(t > 0) only for a restricted class of potentials. It seems 
that in general (even for fast decaying potentials) the 
QNMs appear later (t > t (x,x') > \x — x'\) than in an 
optimal case (t > to(x,x') = \x — x'\). 

There remain following open problems. First, the 
analytic structure (singularities, cuts, asymptotics) of 
G(s, x, x') in the complex plane of s and its dependence 
on the properties of the potential V(x) (asymptotics, an- 
alyticity) need further study. Second, the earliest time 
to(x,x') of validitiy of the QNMs expansion and its de- 
pendence on the shape of the potential is in general un- 
clear. They are connected through the asymptotic be- 
haviour of G(s,x,x') for big s but the relations seem 
very complicated. 

An interesting fact, worth further study, is the appear- 
ance of an additional series of modes in the intermediate 
region before the QNMs appear. However, it has been 
observed only in our example and it is not clear at all if 
this phenomenon is common. 



Acknowledgments 

I want to express my gratitude to Piotr Bizon for in- 
troducing me to the issue of quasinormal modes. 



APPENDIX A: NORMAL MODES, V = 

The homogeneous ODE © has the following solutions 



4>-(s,x) 
(j)+{s,x) 



e sx — e 2sa e ~ sx 



e sx _ e 2sb e -sx 



(Al) 

(A2) 



satisfying the boundary conditions 

<£_(s,a) = 0, (f>+(s,b) = 0. (A3) 
Their Wronskian equals 

W[^J + ]( S ) = ~2 S [e 2sa ~e 2sb ]. (A4) 
The Laplace transformed Green's function reads 



G(s, X, x 1 ) = 



2 s [ e s(a-f)) _ e s(6-o)] 



(A5) 

Its poles, the zeros of the Wronskian W(s), are placed at 



niri 

s n = -, , for 



b — a 

The explicit form of the modes 



n G Z. 



4> n (x) = 0_ (s n , x) = cj> + (s„, x) 



(x) = exp 



sin 



(A6) 



(A7) 



(A8) 



The time-dependent Green's function, obtained by the 
inverse Laplace transform, takes the explicit form 



G(t, x, x') — 

sin ( j^mr) ■ sin ( 4-^-1 

HIT 



° XP 1 



(A9) 



As we know from the alternative way of calculating the 
inverse Laplace transform, the series converges to zero for 
all t < \x — x'\. The series can be formally summed up 
and expressed as a sum of Heaviside theta functions of the 
parameters t, x, x' and hence G(t, x, x') as a sum of Dirac 
delta functions - the result of free propagation inside the 
interval (a, b) with multi-reflections at the boundaries. 



APPENDIX B: NORMAL MODES, V ^ 

The homogeneous ODE, analogous to ©, but with 
non- vanishing V(x), has the following approximate solu- 
tions 



(j>-{s,x) 



e s{x-a) , a __ 



4>+{s,x) 

^s(x-b) 



(l + Ois- 1 )) 



(l + Ois- 1 )) 

(Bl) 
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satisfying the Laplace transformed boundary conditions Their Wronskian equals 



/3 4>-(s,a) - a -g— (s,a) = 

6$ + {8,b)--Y^±-(a,b) = 0. 
with a, f3, 7, 5 € K.. Their Wronskian is 



W[<f>-,4>+](8) 



e s(b-a) , a 



si \ s 



( B2 ) ^[0_^ + ](,s)- S [e si (7 + l) + e- si ( 7 -l)]. (C4) 

(B3) 

The Laplace transformed Green's function reads 
G(s, x, x') — 

[e sx < - e~ sx <] [ e < x >~ L ) ( 7 - 1) + e -^>-L) ( 7 + i)] 



>*<■ (a -jj I 7 + s 



(l + OGs- 1 )) (B4) 



The Laplace transformed Green's function reads 



2s [e sL (j + 1) + e- sL (7 - 1)] 

(C5) 

Its poles, the zeros of the Wronskian W(s), are placed at 



G(s,x,x') = — 
2s 



s(x < -a) I . P \ , -s(K<-a) ( ^ 
e s{*>-b)( 1+ 5 _\ +e -s{*>-b) ( 1 _ 5 _ 



e s{b - a) | a 



(l + OGs- 1 )) (B5) 



sn = ^(^) + n 4, for n.Z. (08) 



The explicit form of the modes 



4> n {x) = 4>^{s n ,x) = 



(-1)" 



^+(s n ,ar) (C7) 



is 



APPENDIX C: LEAKY CAVITY 

For the leaky cavity the homogeneous ODE ijjjj has the 
following solutions 



4>n{x) = 2 sinh(s„a;) = 

/ • \ n ( X / L ) 

( nmx \ / 1 — 7 



1 + 7 

1-7 



(x/L) 



(C8) 



0+(s,x) 



.(s,a;) = e s 



(CI) 



( 7 - 1) + e- s ^- L ) (7 + I) 



satisfying the boundary conditions 

$-(8,0) = 0, (C2) 
s $ + (s,L)+j^±(s,L) = (03) 



The time-dependent Green's function, obtained by the 
inverse Laplace transform, takes the form 



sinh(s„a;) ■ sinh(s n x') 



exp(s„t). (C9) 



As we again know from an alternative way of calculat- 
ing the inverse Laplace transform, the series converges to 
zero for alH < \x — x'\. 
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